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ABSTRACT 

We identify three Kepler transiting planets, Kepler-7b, Kepler-12b, and Kepler-41b, whose orbital 
phase-folded light curves are dominated by planetary atmospheric processes including thermal emis¬ 
sion and reflected light, while the impact of non-atmospheric (i.e. gravitational) processes, including 
beaming (Doppler boosting) and tidal ellipsoidal distortion, is negligible. Therefore, those systems 
allow a direct view of their atmospheres without being hampered by the approximations used in the 
inclusion of both atmospheric and non-atmospheric processes when modeling the phase curve shape. 
We present here the analysis of Kepler-12b and Kepler-41b atmosphere based on their Kepler phase 
curve, while the analysis of Kepler-7b was already presented elsewhere. The model we used efficiently 
computes reflection and thermal emission contributions to the phase curve, including inhomogeneous 
atmospheric reflection due to longitudinally varying cloud coverage. We confirm Kepler-12b and 
Kepler-41b show a westward phase shift between the brightest region on the planetary surface and the 
substellar point, similar to Kepler-7b. We find that reflective clouds located on the west side of the 
substellar point can explain the phase shift. The existence of inhomogeneous atmospheric reflection in 
all three of our targets, selected due to their atmosphere-dominated Kepler phase curve, suggests this 
phenomenon is common. Therefore it is likely to be present also in planetary phase curves that do 
not allow a direct view of the planetary atmosphere as they contain additional orbital processes. We 
discuss the implications of a bright-spot shift on the analysis of phase curves where both atmospheric 
and gravitational processes appear, including the mass discrepancy seen in some cases between the 
companion’s mass derived from the beaming and ellipsoidal photometric amplitudes. Finally, we dis¬ 
cuss the potential detection of non-transiting but otherwise similar planets, whose mass is too small to 
show a gravitational photometric signal but their atmosphere is reflective enough to show detectable 
phase modulations. 

Subject headings: planetary systems — stars: individual (Kepler-7, Kepler-12, Kepler-41) — tech¬ 
niques: photometric 


1. INTRODUCTION 

The availability of space-based high-quality time se¬ 
ries photometry over the last several years has allowed 
the monitoring of transiting planetary systems not only 
during transit or occultation (secondary eclipse) but 
throughout their entire orbit, for a growing sample of 
plan ets. This was done at first by Sy i tzer in the infrared 

('e.g..lKnutson et'ani2007ll2009l[20ialCowan et al.l[2012t 

iLewis et al.ll2013ir ~thus facilitating an improved under¬ 

standing of processes taking place in planetary atmo¬ 
spheres, such as winds that shift the hottest region on 
the planet surface eastwa rd from the substellar point 
(|Showman fc Cuillot] 120021 1. Later, the CoRoT and Ke¬ 
pler space missions enabled similar monitoring in visible 
light (optical), which is sensitive also to light from the 
host star reflected by the planetary atmosphere. 

Moreover, such high-quality space-based optical 
data is sensitive also to non-atmospheric processes 
induce d by the planet -s tar g ra vitational inte r action 
le.g.. iMazehfc Faig ^ 120101 : LShuorer et al.l 120111 : 
iFaigler &: Mazehll2011 L including beaming (aka Doppler 
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boost i ng; e.g., iLoeb fc Gaudil l2003t iZucker et ahl 
2007t iS hporer et M l201d : Ivan Kerkwiik et al.l l201(Tl 
Bloemen et al.l 12011 ) and tidal ellipsoidal distortion 
modulations (w hich we refe r to hereafter simply a s 
ellipsoidal; e.g.. iMor^ 119851 ; iMorris fc: Naftilaiil Il993ll . 
Since both processes are gravitational in origin they 
are both sensitive to the orbiting companion’s mass. 
Therefore, phase curves in the optical allow the study 
of both the planetary atmosphere and the gravitational 
i nteraction between the pl anet and the host star. 

IFaigler fc Mazehl (1201111 proposed to utilize phase 
curve modulations of non-transiting systems to infer the 
existence of an unseen orbiting companion (substellar or 
stellar), by simultaneously modeling all processes, atmo¬ 
spheric and gravitational, and using the latter to esti¬ 
mate the companion’s mass. Their simplistic model con¬ 
sists of simple approximations for each proc ess, and their 
approach h as alrea dy proven successf ul (|Faigler et al.1 
120121 [201^ see also iSIiporer et al1l20lHl . 

Recently, the number of transiting planets with a mea- 
sured optical phase curve has grown significa ntly (e.g., 
lEsteves et al.l 120131 ; lAngerhausen et ahl 1201411 , and the 
existence of a shift between the substellar point and the 
brightest regio n on the planet surface was suggested for 
a few of them (iDemorv et al.ll2013t lEsteves et al.11201^ 
IFaigler fc Mazehl 12015(1 . Such shifts include both east¬ 
ward and westward bright-spot shifts, corresponding to 
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a pre-occultation and a post-occultation maximum (re¬ 
spectively) in the phase curve’s atmospheric component. 

However, analysis of the increasing sample of phase 
curves has resulted in a few inconsistencies. Most no¬ 
tably is a discrepancy that appears in some cases between 
the companion’s mass derived from the beaming photo¬ 
metric amplitude and the mass derived from the ellip¬ 
soidal photometric amplitude. For systems where radial 
velocities (RVs) of the host star are available the RV mass 
confirms the beaming mass in some cases while in others 
it confirms the ellipsoidal mass. This mass discrepancy 
appears in both sta r -planet systems, e. g ., Kepler-13A 
(iShpqrer et al.l I201H iMazeh et all l2012t festeves et al.l 


l2niilHShDorer et al.ll2ni4n. TrES-2 (IBarclav et al.llmli 

and HAT-P-7 (|Esteves et al.l 1201311 . and stellar binaries 
wit h a transiting white dwarf companion, including KOI- 
74 (jvan Kerkwiik et al.ll201C : iBloemen et al.l|201^ . and 


KIC 10657664 |Carter et al.lf2nilh . This points to a gap 
in the understanding of one or more of the processes gov¬ 
erning optical phase curves. 

Confronting the mass discrepancy described above can 
be done by studying phase curves showing only one of the 
processes govern ing their shape. For the beaming effect 
this was done bv iShporer et al.l (j201fl[ l who measured the 
optical (SDSS -5 band) phase curv e of a detached eclips- 
ing double white dwarf binary (jSteinfadt et al.l 1201011 
where the objects’ compact nature results in only the 
beaming effect being detectable at a high signal to 
noise ratio. Here we carry out a study of planetary 
atmospheres using Kepler phase curves that are com¬ 
pletely dominated by atmospheric processes, including 
both thermal emission and reflected light, while the grav¬ 
itational processes (beaming and ellipsoidal) are much 
smaller in photometric amplitude, by an order of magni¬ 
tude or more, and are close to the noise level. Therefore, 
such phase curves present an opportunity for studying 
planetary atmospheres as they are clean in the sense 
that they allow a careful study of the planetary atmo¬ 
sphere without being hampered by the approximations 
done in the simultaneous m odeling of all processes a ffect- 
ing the phase curve shape (jFaigler fc Mazed 1201 111 . For 
the atmosphere-do minated phase cu rves we identify we 
use the formalism of lHu et all (|2015ll for a more detailed 
study of the planetary atmospheres than done previously. 

Out of the ^20 currently known hot Jupiters with 
a measured mass in the Kepler field we have identi¬ 
fied three systems whose phase cu rve is dominated by 
atmospheric process es: Kepler-7 (iLatham et al.l l201fll : 
iDemorv et al.l 1201111. Kepler-12 (iFortnev et al.l 1201111. 
and Kepler-41 ( Santerne et al.l 120111 : iQuintana et al.l 
1201311 . For those three systems the gravitational pro¬ 
cesses are at the 1 part per million (ppm) level in rel¬ 
ative flux, which is more than an order of magnitude 
smaller than the observed phase curve amplitude and 
well within the phase curve noise level. Another five 
systems include significant non-atmospheric processes in 
their phase curve (TrES-2, HAT-P-7, Kepler-13, Kepler- 
76, Kepler-412). For the rest a phase curve signal was 
not detected, either due to insufficient signal to noise ra¬ 
tio, or due to stellar activity, from the planet host star or 
another star within the Kepler pixel mask, that prevents 
identifying the orbital signal. 

We describe our Kepler data analysis and model fit¬ 
ting in Section 2 and show our results in Section 3. In 


Section 4 we discuss our findings and their implications, 
and bring a short summary in Section 5. 

2. DATA ANALYSIS 
2.1. Data preprocessing 

We use here Kepler long cadence data from quarters 2 
through 170. Quarters 0 and 1 are ignored since they 
show a larger scatter than the other quarters. We have 
ap plied several preprocessing steps to the data, follow¬ 
ing [ShporeF^F^ (I2011L1201^ . to prepare it for the at¬ 
mospheric model fitting. We first remove instrumental 
signals, or trends, by fitting the first four cotrending ba¬ 
sis vectors (CBVs) to the data of each quarter using 
the Pyke Python package (|Still &: Barclavl [201211 . We 
checked that increasing to six CB Vs gave consistent re- 
sults well within la, consistent with lShnorer et ^ (I2011L 
l2(rn and IHn et a, 1.1 (jM!). 

We then fit each continuous data segment with a 
5th degree polynomial while ignoring in-transit and in- 
occultation data, and then dividing all data within the 
segment by that polynomial. This step does not affect 
the sinusoidal modulation along the orbit since the du¬ 
ration of each continuous segment is at least an order of 
magnitude longer than the orbital period. Fitting was 
done while iteratively rejecting 5a outliers until none are 
identi fied using the formalism described in lShporer et al.l 
(12014 see their Section 2). Using polynomial of degrees 4 
and 6 did not change the results. Using smaller degrees 
showed an increased overall scatter in the light curves 
while showing the same sinusoidal amplitude, where the 
latter is the astrophysical signal we seek to characterize 
here. Using larger degrees decreased the sinusoidal am¬ 
plitude resulting from the excessively aggressive detrend¬ 
ing. The phase curves before and after the polynomial 
detrending are presented in Appendix 

Next we phase folded each light curve using the known 
ephemeris and binned it using 200 phase bins. We re¬ 
jected 4tT outliers from each bin, which contained about 
230 and 300 individual Kepler measurements for Kepler- 
12 and Kepler-41, respectively. For each bin we calcu¬ 
lated the median relative flux value, and the median 
absolute deviation (also known as MAD) times 1.4826 
to be the uncertainty. Therefore the bins’ uncertainties 
are independent of the error bars assigned to the Kepler 
individual measurements and is based on the true scat¬ 
ter in each phase bin. We compared the measured scat¬ 
ter within each phase bin to the expected scatter based 
on the error bars of the individual measurements and 
Poisson statistics. We found that on average the mea¬ 
sured scatter is 11% and 49% larger than the expected 
scatter for Kepler-12 and Kepler-41, respectively. For 
Kepler-41 the increased measured scatter can be at least 
partially attributed to stellar activity (see Section 4.3). 
Although, for Kepler-12 there is no evidence for stellar 
activity (see Section 4.3) so the increased measured scat¬ 
ter is attributed to residual correlated noise, which in 
principle may exist also in the Kepler-41 data but for 
the latter the stellar activity signal dominates. 

Next we checked whether combining the data from 
different Kepler quarters follows Poisson statistics or 

® Kepler-12 was located on Kepler CCD module 3 during quar¬ 
ter 4, when that module failed. So it is missing part of quarter 4 
and the entire quarters 8, 12, and 16. 
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whether it adds some noise due to, e.g., instrumental dif¬ 
ferences between quarters or data detrending. For each 
object we generated a phase folded and binned light curve 
for each quarter. We then calculated for each phase 
bin the ratio between the average uncertainty, averag¬ 
ing across all quarters, and the uncertainty of the phase 
bin when using all Kepler data, and further dividing by 
the square root of the number of quarters. For both 
the Kepler-12 and Kepler-41 data sets the average ratio 
across all phase bins was within 1-2 % from unity, show¬ 
ing that combining data from different quarters closely 
follows Poisson statistics. 

Finally, we checked that our results do not depend on 
the phase bin width by analyzing also phased light curves 
with 100 and 400 bins and confirming the results are 
identical. 

Before moving to the atmospheric modeling we re¬ 
moved from the phase-folded light curves the orbital 
modulations induced by beaming and el lipsoidal, based 
on the known planet mass in each system (jSanterne et al.l 
l2011l : lDemorv et aI1l2011h . This is done for completeness 
since as mentioned above, for the systems analyzed here 
these modulations are small, at the noise level, and the 
phase curve amplitude is dominated by atmospheric pro¬ 
cesses (see Figure [TJ so any approximations in determin¬ 
ing the shape of the beaming and ellipsoidal phase curves 
do not affect the overall phase curve shape significantly. 

2.2. Phase curve analysis with atmospheric modeling 

We analyze the Kepler phase curves using the inter¬ 
pretation framework described in I Hu et al.l ( 2015f) . The 
model efficiently computes reflection and thermal emis¬ 
sion contributions to the phase curve, considering both 
hot spot shift due to equatorial winds and inhomoge¬ 
neous atmospheric reflection due to patchy clouds. The 
model assumes the atmosphere has clear longitudes and 
cloudy longitudes, and the cloudy part is more reflective 
than the clear part. The model assumes the cloud dis¬ 
tribution is controlled by the temperature distribution, 
which is in t urn approxima t ed by the analytical model 
proposed bv ICowan fc Agoll (|201ir) . The model also as¬ 
sumes the thermal emission comes from a thermal photo¬ 
sphere that may be hotter than the equilibrium tempera¬ 
ture. Here, the equilibrium temperature, as a function of 
longitude, is defined to be the atmospheric temperature 
without any greenhouse effects, and is calculated from 
the Bond albedo and the heat redistribution efficiency. 

We fit the phase curve to this model characterized by 
hve parameters: the Bond albedo, a heat redistribution 
efficiency, a greenhouse factor, the condensation tem¬ 
perature of clouds, and a reflectivity boosting factor by 
clouds. The Bond albedo describes the overall reflectiv¬ 
ity of the planet, and takes any value between 0 and 1. 
The heat redistribution efficiency, the ratio between the 
radiative timescale and the advective timescale, controls 
the longitudinal temperature distribution. The magni¬ 
tude of this parameter describes how well heat is trans¬ 
ported: the larger the parameter, the smaller the longi¬ 
tudinal temperature variation becomes. The sign of this 
parameter describes the direction of the equatorial winds: 
a positive sign indicates super-rotating eastward winds, 
and a negative sign indicates westward winds. The green¬ 
house factor is the ratio between the temperature of the 
thermal photosphere and the equilibrium temperature. 


and takes any value not smaller than 1. The cloud con¬ 
densation temperature, another free parameter, deter¬ 
mines the longitudinal boundaries of the clouds. Finally, 
the reflectivity boosting factor by clouds is the propor¬ 
tional increase in reflectivity of an atmospheric patch 
when it becomes cloudy, and takes any positive value. 
The rest of the syst em parameters are known and taken 
from the literature (Ifirtnev et al.l 120111: iSanterne et al.l 
120111 : iQuintana et al.ll2013lb 

We use the Markov-Chain Monte Carlo (MCMC) 
method (|Haario e t al.| |2006l) to explore t he parameter 
space. Following IGelman fc Rubinl (|1992D we calculate 
two Markov chains for each of the two objects, each 
chain containing 10® steps where the first half of each 
chain is considered the “burn-in” period and is removed 
once the chain is completed. We verify that the R values 
for all parameters are less tha n 1.01 to ensure conver¬ 
gence (|Gelman fc Rubinl 119921) . and that in both cases 
the results from the two chains are identical. The al¬ 
lowed ranges and the prior distributions of the model 
parameters are as follows: the Bond albedo uniformly 
ranges in [0,1]; the heat redistribution efficiency uni¬ 
formly ranges in [0,100]; the greenhouse factor uniformly 
ranges in [1,2]; the cloud condensation temperature uni¬ 
formly ranges in [1000, 3000] K; and the reflection boost¬ 
ing factor uniformly ranges in [0,100]. The redistribu¬ 
tion efficiency parameter is only allowed to take a pos¬ 
itive value because atmospheric cir culation models pre¬ 
dict eastward equatorial winds (e.g. JShowman &: Guillod 
120021 : IShowman &: Kasml I2013H . This choice effectively 
requires any post-occultation phase offset to be explained 
by the asymmetric reflection components induced by the 
patch y clouds, rather than the thermal emission compo¬ 
nent (|Hu et al1l2015fl . Based on the converged Markov 
chain, we derive the occultation depth, the phase curve 
amplitude, and the phase shift of the phase curve maxi¬ 
mum from occultation phase. The latter is defined to be 
positive for post-occultation maximum. 

After carrying out the analysis once, we noticed that 
for Kepler-12b the reduced of the residuals was 1.12. 
Therefore we repeated the analysis of that phase curve 
while increasing the bins’ uncertainties to bring the re¬ 
duced to unity. This has increased the uncertainty 
on the htted parameters by up to a few percent. For 
Kepler-41b the reduced x^ was 0.86 in the first analysis, 
so no further analysis was done. 

To test our results we applied the so-called “prayer 
bead” approach to both objects, where the residuals are 
cyclicly permuted and added back to the model, and the 
new data set is then refitted. This approach preserves 
any correlated noise features in the data and produces 
a distribution for each fitted para meter from which its 
uncertainty can be estimated (e.g., IBouchv et ahl 120051 : 
ISouthworthI 120081 : IWinn et al.l 1200^ . For both objects 
the prayer bead analysis did not result in larger uncer¬ 
tainties than the original MGMG analysis, so we use the 
latter as our final result. 


3. RESULTS 


Our resulting fitted light curve models are shown in 
Figure [T] and the fitted and derived parameters are listed 
in Table [II In the latter w e list also the results for Kepler- 
7b from Hu et all ([2015f) for comparison and complete¬ 
ness. In Appendix iBl we show the correlations between 















































4 


Shporer & Hu 


the five fitted model parameters for each of the two ob¬ 
jects. 

For both Kepler-12b and Kepler-41b we find that their 
phase curves are dominated by atmospheric reflection 
and we identify a statistically significant shift of the 
phase curve maximum from the occultation phase where 
the maximum is post-occultation, similar to Kepler- 7 b 

(|Demorv et al.l201^lHu et al.ll201alWebber et al.ll2015ll . 


We fitted Kepler-41b with the same simple paramet¬ 
ric model as we did for Kepler-12b, consisting of a si¬ 
nusoidal amplitude and phase. We derived a peak to 
peak sinusoidal amplitude of 39.2±4.2 ppm, and a phase 
of maximum shifted later than the occultation phase by 
25.2 ± 5.9 deg. These are within 1.5 (t from the results of 
our original analysis and with similar uncertainties (see 
Table P. 


3.1. Kepler-12b 

We measure a phase curve amplitude of 17.5± 2.2 ppm 
(part per million) and an eclipse dept h of 14.3 ± 2.5 ppm, 
consistent with previou s studies (|Angerhausen et ^ 
I20l4lfeteves et al]l2015ll but with smaller uncertainties. 
We detect a clear phase curve asymmetry, with a phase 
curve maximum shifted by g deg from occultation 

ph ase. This is consiste nt with the phase shift reported 
by lEstev es et al .l (l2015f) b ut Ri 2tT smaller than that of 
lAngerhausen et al.l (|20l4 D. Angerhausen, private com¬ 
munication). We find that the phase curve amplitude 
is greater than the eclips e depth by ^Igi in agreement 
with the measurements of lEsteves et al.l (|2015[ 1. The dif¬ 
ference between the phase curve amplitude and eclipse 
depth is enhanced by the relatively large phase offset. 

To further test our analysis we carried out a separate 
analysis where we fitted a simple parametric model in¬ 
cluding a sinusoidal amplitude and phase. This is similar 
to previous a nalyses (see a summ ary of previous results 
in Table 9 of lEsteves et al.l[2015H . We derived a peak to 
peak amplitude of 19.0 ± 2.2 ppm and a phase of maxi¬ 
mum shifted later than the occultation phase by 58.5±6.8 
deg. These values are ^ la away from that of our original 
analysis and with similar uncertainties (see Tabled I. 

Looking at Kepler-12b phase curve (Figure [T] left 
panel), there are a few consecutive bins around phase 
0.2 that fall below the fitted model. While it could be an 
astrophysical signal that is not included in our model, it 
could also be a statistical, correlated noise feature. Al¬ 
though, correlated noise is expected to average out in the 
pha se fo lded light curve. Either way, as described in Sec¬ 
tion [22] we have inflated the bins’ uncertainties to bring 
the residuals reduced to unity, and the prayer bead 
analysis did not result in larger error bars for the fitted 
parameters compared to the MCMC analysis. 


3.2. Kepler-4 lb 

Although an asym metry in Kepler-41 b phas e curve 
is not identi fied bv lAngerhausen et al.l (j2014ll . it is 
identified bv lEsteves et al.l (1201^ and is visually ap- 
parent in iSanterne et aLl (|2011l . their Fig. 3) and 
iQuintana et al.l ( 20131 their Fig. 2). Here we confirm 
the post-occultation phase maximum which we mea¬ 
sure to be shifted by 17.8 ± 6.3 deg. We derive an 
eclipse depth of 49 . 1 I 5 2 ppm and a phase amplitude 
of 48.7 ± 5.2 ppm, within l-2a from previous results 
but w ith typically smaller e rror bars (ISanterne et al.l 
2014 [Quintana et al.l 120 131: lAn gerhausen et al.1 12014 
Esteves et al.ll2015[ l. lEsteves et al.l ()2015l l report a phase 
amplitude greater than the eclipse depth (although at 
low statistical significance of only 1.2cr), while here we 
find a phase curve amplitude that is highly consistent 
(within O.lcr) with the eclipse depth, indicating either 
poor heat redistribution or large planetary albedo. 


4. DISCUSSION 
4.1. Atmospheric constraints 

Our analysis of Kepler-12b and Kepler-41b Kepler 
phased light curves confirms and refines previous results 
and includes a more detailed atmospheric analysis than 
done previously. 

We identify an asymmetry in the Kepler phased light 
curves for the two planetary systems studied here, where 
the maximum is post-occultation. Our interpretation is 
that this asymmetry is due to the brightest region on 
the planetary surface being shifted westward from the 
substellar point, that in turn originates from an asym¬ 
metric surface distribution of clouds in the planetary at¬ 
mosphere. The post-occultation maximum requires that 
the reflective clouds exist on the west side of the substel¬ 
lar point but not on the east side. Such asymmetric cloud 
distribution may be a result of eastward heat transport 
driven by zonal jets in the atmosphere. In addition, the 
clear part of the atmosphere must be quite dark, prob¬ 
ably due to absorption by Alkali metal vapors deep in 
the atmosphere, whereas th e cloudy part appears to b e 
bright due to Mie scattering (iSudarskv et al.t2000ll2003ll . 

For Kepler-41b, in order to produce a phase varia¬ 
tion consistent with Kepler observations, the east side 
of the sub-stellar meridian should be devoid of reflec¬ 
tive clouds and have a reflectivity less than 0.2 (Icr limit, 
see Table [H). Figure [ 1 ] shows that the thermal emission 
component has the opposite phase offset than the asym¬ 
metric reflection component. On the west side, where 
reflective clouds are present, the reflectivity should be 
greater than 0.3. The overall Bond albedo of the planet 
is constrained to be 0.25 ± 0.06. The cloud condensation 
temperature parameter is found to be 1580^120 
compare with material properties, one should also con¬ 
sider that clouds can form dee p in the atmospher e and 
hence at a higher temperature (|Hu fc Seageill20l4l . Af¬ 
ter multiplying by the derived greenhouse factor the Icr 
cloud condensation temperature range is 1530-1910 K. 
Silicates (MgSiOa and Mg 2 Si 04 ) are consistent with this 
temperature range and are highly reflective, so are there¬ 
fore candidate condensates. 

Comparing Kepler-12b with Kepler-41b, we find that 
Kepler-12b is a more extreme case. Its phase curve 
is completely dominated by asymmetric surface distri¬ 
bution of reflective clouds, whereas the phase curve 
of Kepler-41b is dominated by symmetric reflection 
(see Figure [T|). Because the cloud distribution is con¬ 
trolled by the cloud condensation temperature, the post- 
occultation phase offset yields a tight constraint on the 
latter (I 5 O 6 I 24 K). Multiplying by the derived green¬ 
house factor the Itr cloud condensation temperature 
range is 1480-1670 K. To produce the phase offset, the 
clear part of Kepler-12b atmosphere needs to have an 
albedo very close to zero, lower than that of Kepler-41b. 
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Figure 1. Phase folded and binned Kepler light curves of Kepler-12b (left) and Kepler-41b (right). Top panels show relative flux (planet 
flux, Ep, relative to the in-occultation stellar flux, Fs) in part per million (ppm) versus orbital phase. Occultation (secondary eclipse) is at 
phase 0.5 and transit (not shown) at 0.0. The black curve shows the best-fit model phase curve and the colored lines show contribution of 
thermal emission (red), symmetric reflection (green), and asymmetric reflection (blue). The model phase curve has a significant contribution 
from the asymmetric reflection component due to patchy clouds. The clouds concentrate on the west side of the substellar point, due to a 
hot spot shifted eastward. The bottom panels show the residuals, data subtracted by the model, with a gray dashed line at zero residual 
for reference. For these systems the orbital gravitational processes, beaming and ellipsoidal, are typically at the 1 ppm level, that is close 
to the noise level and much smaller than the overall phase curve amplitude which is dominated by atmospheric processes. 


This can be naturally explained as a higher metallic- 
ity and therefore greater Na and K absorption in th e 
atmosphere of Kepler-12b (jSudarskv et alJl2000l I2003D . 
However, the host star in the Kepler-12 system is rela¬ 
tively metal-poor compare d to Kepler-41 (jFortnev et alJ 
1201 It iSanterne et alJl2nilh , therefore it is unclear why 
Kepler-12b would have a more metal-rich atmosphere 
than Kepler-41b. 

Overall, our analysis shows that it is reasonable to as¬ 
sume that an asymmetric cloud distribution is a common 
phenomena for hot Jupiters, but the details for the clear 
and cloudy parts of their atmospheres can still manifest 
great diversities. 

For both planetary atmospheres the current data con¬ 
strain the Bond albedo and the cloud condensation tem¬ 
perature, and put an upper limit on the greenhouse fac¬ 
tor. Among these parameters, we see that the Bond 
albedo and the cloud condensation temperat ure are 
somewhat correlated, as expected (See Figures IBII and 


IB2|) . This is because to produce a patchy cloud solution 
the cloud condensation temperature needs to be between 
the maximum and the minimum temperature of the day- 
side, both affected by the Bond al bedo. The sam e cor¬ 
relation is observed for Kepler-7b (jHu et al.ll2015H . We 
note that the efficiency of heat redistribution cannot be 
sufficiently constrained by the current data. Although 
the best-fit models suggest efficient heat redistribution, 
many more models would be permitted at an expense of 
minimally increasing For example, we find that much 
less efficient heat redistribution and a much greater re¬ 
flectivity contrast between the cloudy part and the clear 
part could ren der an almost-as -good fit for Kepler-41b. 
As discussed in iHu et al.l (|2015f) . this insensitivity to the 
exact value of the heat redistribution efficiency is tied to 
the uncertainty in the exact condensation temperature 
of the cloud material. 

The atmospheric analysis done here for Kepler-12b and 
Kepler-41b and the analysis done for Kepler-7b by others 
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Table 1 


Parameter 

Kepler-12b 

Kepler-41b 

Kepler-Tb® 

Fitted parameters: 

Bond albedo 

0.12 ±0.02 

0.25 ± 0.06 

0.42 ± 0.01 

Heat redistribution 

46/^27 

Q q+36.1 
^•^-1.3 

40 + 33 
^^-29 

Greenhouse factor 

1.05 ±0.05 

l-O+aJi 

1 08^^°^ 
J--^c>_0.05 

Cloud condensation temp. [K] 

I 5 O 6/24 

1580/(50 

1480 ± 10 

Cloud refiectivity factor 

45+36 
^0 27 


28//( 

Derived Parameters: 

Eclipse depth [ppm] 

14.3 ± 2.5 

49 . 1 / 5.2 

39.0/(J 

Phase amplitude [ppm] 

17.5 ±2.2 

48.7 ±5.2 

46.6 ± 1.3 

Phase offset [deg] ^ 

A7.6t.fl 

17.8 ±6.3 

36.2 ± 2.7 

West longitude [deg] ^ 

-33± 15 

09 1 +cJi.c$ 
OZ.i_i 2 3 

-11.2 ±2.7 

East longitude [deg] ^ 

90 

90 

90 

A , 

■^cLear 

u.uuyi_Q QQ 5 ^ 

0.026/0;J(4 

U.UO+_Q Q20 

4 , ^5 

■^cLoua 

'-'•^^-0.14 

0 . 68 /°j( 

0.92 ± 0.04 


Note. — 

1: Defined to be positive for post-occultation maximum. 

2: Longitudinal boundary of clouds in the western hemisphere, with the 
substellar point defined as zero longitude. 

3: Longitudinal boundary of clouds in the eastern hemisphere. 90 degree 
means the eastern hemisphere is completely cloud-free. 

4: Refiectivity of the clear longitudes. 

5: Refiectivity of the cloudy lon gitudes. 

6 : Results from [Hu et al.l l|2015l'l . 


(|Demorv et al.l[2Ml : iHu et al.]l2015l : iWebber et al.ll2015[) 

point to a similar scenario, that the optical Kepler phase 
curve is dominated by reflected light (as opposed to ther¬ 
mal emission) and that the shift between the phase curve 
maximum and planetary occultation phase is a result of 
a non-uniform longitudinal cloud coverage. The irradia¬ 
tion received by these three planets is similar, and hence 
the cloud condensation temperatures derived from our 
phase curve analysis are also similar. To form the patchy 
clouds prominent longitudinal distribution the conden¬ 
sation/evaporation timescale of the cloud parti cles needs 
to be shorter than their advection timescale (|Hu et al.l 
1201511 . This condition is uniquely met by hot Jupiters 
that have fast zona l winds in their atmospheres (e.g. 
iKnutson et al.l[2?i(y^ . 

4.2. Phase shift implications 

The three planets identified here as having a phase 
curve dominated by atmospheric processes, hence they 
allow a direct view of their atmospheres, all show a 
phase shift in the form of a post-occultation maximum. 
The identihcation of phase shifts in these atmosphere- 
dominated phase curves suggests they are common and 
exist also in phase curves of other hot Jupiters where they 
are more difficult to identify uniquely as those systems 
show additional phase modulations, induced by beaming 
and tidal ellipsoidal distortion. This has several potential 
implications we discuss below. 

One obvious implication of the phase shift is on atmo¬ 
spheric studies using only the occultation and not the 
entire phase curve. Such studies commonly assume that 
the occultation depth is a flux measurement of the bright¬ 
est region on the planetary surface, which is not the case 
for shifted phase curves. In those cases, the secondary 
eclipse depth underestimates the planetary brightest- 
region flux, which in turn can lead to underestimating 
the planetary reflectivity and/or thermal emission. 

Another implication of the phase shift involves the 


BEER model (|Faigler &: Mazehl 1201 ID . developed to 
search for non-transiting systems through the detection 
of photometric orbital phase modulations. It assumes 
that all processes, including gravitational (beaming and 
ellipsoidal) and atmospheric (reflected light and thermal 
emission), are all well aligned with the phase of transit 
and occultation and that the orbit is circular. Mean¬ 
ing, the extrema of each effect occurs exactly at either 
transit, occultation, or quadrature phases. When the 
reflection component is shifted, with either a post- or 
pre-occultation maximum, then the BEER model will 
give a biased estimate of the beaming amplitude. This is 
so since the beaming and reflection components are the 
sine and c osine components of the same frequency modu¬ 
lation Isee lFaigler fc Mazehll20lTI for more details). This 
in turn will make the companion mass estimate based 
on the beaming amplitude, m 2 ,beam, disagree with the 
estimate based on the ellipsoidal amplitude, m 2 ,eiiip- A 
post-occultation maximum, as identified here, will de¬ 
crease the beaming amplitude measured by BEER, lead¬ 
ing to m 2 ,beam < m 2 ,eiiip, while a pre-occultation maxi¬ 
mum will lead to m 2 ,beam > m 2 ,elUp- 

As noted earlier, the mass discrepancy between the 
beaming-based mass estimate and the ellipsoidal-based 
mass estimate has already been measured for several 
systems, including transiting planets and stellar bina¬ 
ries with a transiting white dwarf companion (see Sec¬ 
tion I). Another source of disagreement between the 
beaming-based and ellipsoidal-based mass estimates are 
the approximations in the interpretation of the ellip¬ 
soida l light curve shape (|MorrisllI985HMorris fc NaftilanI 
Il993t iKahl et al.l 1200811 . especially in cases where the 
primary star rotation is not synchronized with the or¬ 
bit, or its spin axis is not aligned with the orbital angu¬ 
lar momentum axis. Other factors that might be fur¬ 
ther contributing to this disc repancy are poor u nder- 
stating of the host star (e.g., iShporer et al.l 1201^ and 
non-zero orbital eccentricity. Therefore, deciphering the 
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origin of the discrepancy requir es a careful case by case 
study. iFaigler fc Mazehl (1201 51) recently suggested that 
the discrepancy for Kepler-13Ab, TrES-2b, HAT-P-7b, 
and Kepler-76b, identified in Kepler data, can be at¬ 
tributed to a pre-occultation phase curve maximum of 
the atmospheric component. They use the difference be¬ 
tween the two mass estimates to derive the angular shift 
of the brightest region on the planetary surface from the 
substellar point and attribute this shift to asymmetric 
thermal emission of the planetary atmosphere. 

Yet another possible implication of the phase shift and 
the planetary asymmetric surface brightness distribution 
is an asymmetry, or a distortion, of the occultation light 
curve ingress and egress. When identified at sufficiently 
high signal to noise ratio it can be used to map the 
planet’s surfac e brightness, a method known as eclips e 
mapping (e.g., IWilliams et al.ll^OOdUde Wit et al.ll2012D . 
When the distortion is not directly identified it can still 
lead to a shift in the measured occultation phase when 
mode led using a symmetr ic occultation light curve model 
('e.g.. iShporer et al1l2014l) . However, the Kepler data of 
the three transiting planets studied here is not sensitive 
to this minute effect. 


Finally, we note that a longitudinally non-uniform 
atmosphere will affect the shape of the polar i zation 


reflected light phase curve (e.g.. ISeager et al. 

[20?)(Tl: 

iBerdvugina et ai.ll2011l:IMadhusudhan & Burrows 

EbI). 


4.3. Can the atmospheric-induced modulations be used 
to detect similar but non-transiting exoplanets? 

The three exoplanets studied here (Kepler-7b, Kepler- 
12b, and Kepler-41b) all have a mass of roughly half 
a Jupiter mass and orbital semi-major axes within 0.1 
AU. Such a planet mass gives a gravitational signal (due 
to the beaming and ellipsoidal effects) at the ^1 ppm 
level, which is not expected to be detectable in Kepler 
data. On the other hand, their Kepler light curves show a 
clear atmospheric photometric orbital modulation which 
if detectable can reveal the existence of a faint low-mass 
companion even when the system does not transit, as 
discussed in more details below. Therefore, such planets 
represent a class of non-transiting planets that although 
are typically beyond th e reach of the BEER approach 
(|Faigier fc Mazehl 120111) . as the latter is designed to si¬ 
multaneously detect all orbital modulations, can poten¬ 
tially still be identified through their atmospheric photo¬ 
metric orbital modulation resulting from their relatively 
large geometric albedo. The latter is similar to the ap- 
proach discussed in detail bv iJenkins &: Dovli (|2003D . 

Since the overall phased light curve shape is similar to 
a pure sine a simple Fourier-like period analysis could 
be used to identify it, at least in principle. This will 
result in a measurement of the orbital period, and the 
photometric modulation amplitude will give a mutual 
constraint on the companion’s radius, geometric albedo, 
and brightness temperature, at the observed wavelength. 
The companion’s (minimum) mass and the true nature 
of the system can then be obtained vi a radial velocity 
follow -up measurements. Inter estingly, iCharninet et ahl 
(j2011f) and iSilvotti et al.l (|2014[) used a similar approach 
to claim the discovery of non-transiting small planets or¬ 
biting the hot B subdwarf stars Kepler-7[Q and Kepler- 

® Also known as KOI-55, KIC 5807616, and KPD 1943+4058. 


42^1, respectively, at short orbital periods. 

To see whether Kepler-7b, Kepler-12b, and Kepler-41b 
could have been identified as candidate hot Jupiter sys¬ 
tems even if they were not transiting we show in Figy 
ure [2] a Lomb-Scargle (L-S; iLombI 119761 : iScargld 1198^ 
periodogram of their Kepler light curves after removing 
the transit and occult a.tion data. We used th e fast L-S 
algorithm described bv I Press fc RvbicO (I1989D as imple¬ 
mented in the Matlab/plomb function. In all three pan¬ 
els the dashed horizontal black line marks the L-S power 
level corresponding to a false alarm probability of 10“®, 
and the frequency component at the orbital frequency is 
marked by a gray arrow. 

Kepler- 41 (Figure [H bottom p anel) is known as an ac¬ 
tive star (iQuintana et al.l 120131) and its periodogram is 
dominated by the rotational modulation at frequencies 
smaller than the orbital frequency. Yet, there is a clear 
signal at the orbital frequency. For Kepler-7 and Kepler- 
12 (Figure [5] top and middle panels) the strongest fre¬ 
quency components are at the orbital frequency. Mean¬ 
ing, they would have been identified as periodic variables 
at the orbital period even with no transits. In fact, in¬ 
cluding the data during conjunctions, which was removed 
here, would make the orbital frequency component even 
stronger, as we have confirmed through simulations. 

However, there are other processes that may result 
in photometric sinusoidal variability, so the simplistic 
approach desc ribed above may suffer from a high false 
positive rate (jJenkins fc Dovld I2003D . One such sce¬ 
nario involves an active star, where the combination 
of the stellar rotation with the non-uniform spot sur¬ 
face distrib ution results in sine-l i ke photometric modula - 
tions fe.g.. iHartman et aI1l2011l : iMcQuillan et al.ll2014D . 
Although, such modulations vary with time in ampli¬ 
tude and phase, as opposed to a steady orbital signal, 
so a more sophist icated approach, e.g., using wavelets 
(|Bravo et al.ll20l3) . can differ between activity-induced 
modulatio ns and orbital modula tions. Moreover, as 
shown by IMcQuillan et al.1 (|2013D hot Jupiters orbiting 
at short periods are typically hosted by stars rotating 
at significantly longer rotation periods, as is the case 
for Kepler-41b. So the orbital and rotational frequency 
components are expected to be well separated in the pe¬ 
riodogram. 

Stellar asteroseismic pulsations might also be a false 
positive scenario, although, typical pulsations of main 
sequence stars are at a much shorter time scale than h ot 
Jupiters orbital periods ('e.g.. iChanlin fc Miglioll20T^ . 

Yet another false positive scenario can originate from 
blending of nearby stars on Kepler pixels. In that sce¬ 
nario the light from a faint variable star, for exam¬ 
ple an ellipsoidal variable (contact stellar binary), can 
be blended with that of a brighter non-variable star to 
give an observed low-amplitude sinusoidal modulation. 
This scenario, which presents a problem for any vari¬ 
ability study using Kepler data, can be ex amined using 
light curves of the individual Kepler pixels (iBrvson et al.1 
120131) combined with high angular resolution imaging of 
the target that maps the stars within the K epler pixel 
mask ('e.g.. lHorch et al.ll20l4 iLaw et al.ll2014l ). 

Although the detection of hot Jupiters through their 
atmospheric modulations alone is expected to be more 

KIC 10001893. 
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prone to false positive scenarios than BEER, it can al¬ 
low access to (non-transiting) exoplanets in a lower mass 
range, down to below a Jupiter mass, where exoplan¬ 
ets are far more common than at a few Jupiter mass or 
above where they become detectable by the B EER ap¬ 
proach (jShporer et al.l[20lll : lEaierler et al.ll2013[ l. 

5. SUMMARY 

We have identified three Kepler transiting planetary 
systems whose light curves are dominated by atmo¬ 
spheric processes, thus allowing a direct view of their 
atmospheres, independent of any approximations done 
in the simultaneous modeling of additional orbital ef¬ 
fects including beaming and ellipsoidal. We presented 
the analysis of two of them here, Kepler-12b and Kepler- 
41b. while the third , Kepler-7b, was analyzed by other s 

(|Demorv et al.l20T^lHu et al.ll2015l : IWebber et ai.ll20T^ . 

Our analysis of Kepler-12b and Kepler-41b phase 
curves confirms th e light curves’ asymmetric nature 
(|Esteves et al.l[2M5ll where the maximum brightness ap¬ 
pears post-occultation. Our atmospheric modeling shows 
that reflective clouds located on the west side of the sub- 
stellar point can best explain these phase shifts. We 
hypothesize that the asymmetric cloud coverage is the 
result of a longitudinal temperature variation driven by 
zonal jets as the dominant form of atmospheric circula¬ 
tion on those planets. 

The fact that all three systems with atmospheric- 
dominated Kepler phase curves show a post-occultation 
shift of the phase curve maximum suggests it is a com¬ 
mon phenomenon, and similar shifts are likely to exist in 
phase curves of other hot Jupiters where they are more 
difficult to identify uniquely as they contain additional 
phase modulations components, induced by beaming and 
tidal ellipsoidal distortion. 

The phase curve ana l ysis u sing the atmosphere model 
framework of lHu et’abl ()2015h applied here to Kepler-12b 
and Kepler-41b shows that key properties of clouds in 
the atmospheres of highly-irradiated exoplanets can be 
derived from their visible-wavelength phase curves. Fur¬ 
ther understanding can be obtained by the simultaneous 
study of phase curves in other wavelength regimes, for 
example by using Spitzer to obtain IR phase curves, or, 
a dedicated missio n like the proposed FINESSE mission 
(|Deroo et al.ll2012]l to measure phase curves throughout 
a wide wavelength range from the optical to the IR. 

The significant atmospheric photometric orbital mod¬ 
ulation of Kepler-7b, Kepler-12b, and Kepler-41b, along 
with the westward bright-spot shift have possible impli¬ 
cations on the search for non-transiting but otherwise 
similar planets using th e BEER approach in Kepler data 
(|Faigler fc Mazehl 1201111 . and in high-quality photomet¬ 
ric data of other curre nt and f uture space mission s like 
K2 (iHowell et al.ll201^. TE SS l)Ricker et al.ll201'3l . and 
PLATO (IRaiier et al.ll2nil . 

Atmospheric processes (thermal emission and reflected 
light) resulting in an asymmetric surface brightness are 
likely to be at play in phase curves showing also the grav¬ 
itational processes (beaming and tidal ellipsoidal distor¬ 
tion). Such asymmetry can explain some of the mass dis¬ 
crepancy cases, where there is a discrepancy between the 
beaming-based mass estimate and the ellipsoidal-based 
mass estimate (see Section 4.2). Therefore it is a step for¬ 
ward in the detailed understanding of visible-light phase 


curves, and using them through the BEER approach for 
looking for non-transiting faint low-mass binary com¬ 
panio ns in short orbital periods (|Faigler fc Mazel3l201lL 
120151 ). However, asymmetric surface brightness can¬ 
not ex 2 iain_all__the_jmss_^cre 2 ancY_cases(e;£^ KOI- 
74, Ivan Kerkwiik et 'aL 20K; IBloemen et al.l 1 20121: KIC 
10657 664. ICarter et allHoTI : see also iRappaport et al.l 
1201511 . Hence, fully understanding visible-light phase 
curves requires further study, especially in cases of hot 
early-type host stars (like KOI-74 and KIC 10657664) 
where the ellipsoidal distortion is expe cted to behave dif - 
ferently than for cooler late-type stars (|Pfahl et al.l2008ll . 

Finally, we have noted that although non-transiting 
but otherwise similar planets to those analyzed here are 
not expected to show a detectable gravitational signal, 
they might still be detectable through their atmospheric 
photometric orbital modulations. 
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41 (bottom) after removing data during transit and occultation. We used the fast L-S algorithm described by IPress Rvbickl l| 198911 
as implemented in the Matlab/plomb function. In each panel the gray arrow points to the frequency component that coincides with the 
orbital frequency, and the dashed horizontal black line marks the L-S power level corresponding to a false alarm probability of 10~®. 
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APPENDIX 

A. PHASE CURVE FIGURES BEFORE AND AFTER POLYNOMIAL DETRENDING 

Figure IXTl presents the two phase curves analyzed here, of Kepler-12 and Kepler-41, while showing them with (red) 
and without (black) the polynomial detrending (see Section [2]). For Kepler-41 it is visually clear that the polynomial 
detrending decreases the overall scatter, although the difference is minor for Kepler-12. 
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Figure Al. Phase curves with (red) and without (black) the polynomial detrending (see Section (JJ, for Kepler-12 (top) and Kepler-41 
(bottom). Similarly to Figure [T] the light curves present the flux from the planet, Fp, relative to the in-occultation flux from the host star, 
Fb, in part per million (ppm) versus orbital phase. Error bars are not shown to avoid cluttering the figure. 


B. CORRELATION MATRIX FIGURES 

We present in Figures |BT] and |B2] the correlation matrix figures for Kepler-12b and Kepler-41b, respectively, where 
the correlations are shown between any pair of the five fitted model parameters. 
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Figure Bl. Correlation matrix for Kepler-12b showing the correlation between the five fitted model parameters (see Section l2.2ll : Bond 
albedo {Ab), heat redistribution efficiency (e), greenhouse factor (/), clouds condensation temperature (Tc), and clouds reflectivity boosting 
factor (n). 
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Figure B2. Same as Figure [BT] for Kepler-41b. 































